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Sulphur abundance in very metal-poor starsQ 
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ABSTRACT 

We have obtained high-resolution high signal-to-noise spectra of the S I near- 
infrared doublet at 8694 A in eight metal-poor stars with metallicities in the range 
—0.6 <[Fe/H]< —3.0. Elemental sulphur abundances were derived for six targets 



^ . and upper limits were set for two of the most metal-poor stars in the sample. The 

t^ ; sulphur doublet at 8694.62 A has been detected in three stars (HD 19445, HD 2665 

O 

C~^ ' and HD2796) with [Fe/H]< —1.9. Our observations combined with those available in 

o ; 

the literature indicate a monotonic increase of the [S/Fe] ratio as [Fe/H] decreases, 

^ , reaching values [S/Fe]~0.7-0.8 below [Fe/H]=— 2. We discuss plausible scenarios for 



interpretation of these results. 
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Introduction 



The standard scenario for the chemical evolution of the Galactic Halo (see e.g. Tinsley 1979) 
is based on the different yields of Type la and Type II Supernovae. It is well established that 
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Type II supernovae mostly produce the so-called a-elements (O, Ca, Mg, S, Si and Ti) while 
Type la supernovae are responsible for the synthesis of the Fe-group (Fe, Co, Ni, Zn). Given a 
different life-times of both types of supernovae (10^ years for Type II and 10^ years for Type la), 
it is expected that there will be overabundance of a elements at [Fe/H]< —1. At first glance, 
this behavior is observed in many metal-poor halo star (McWilliam 1997). The precise trends 
of the [a/Fe] ratio versus [Fe/H] are still not clearly established. For instance, low [a/Fe] ratios 
were found in several metal-poor halo stars (Carney et al. 1997) and a rather complicated picture 
emerged from the comprehensive work by Idiart & Thevenin (2000) who studied [Ca/Fe] and 
[Mg/Fe] ratios in 252 halo and disk stars. On the other hand, the [0/Fe] ratio in unevolved stars 
show a linear rise (Israelian, Garcia Lopez & Rebolo 1998, Israelian et al. 2001) with a slope 
— 0.33ib0.02 in the metallicity range 0<[Fe/H]< —3. Apparently there is no general agreement 
about a unique plateau [a/Fe]~ const at [Fe/II]< —1 for a-elements such O, Ca and Mg. This 
statement can be extended to another a-element: sulphur. 

In his extensive analysis Frangois (1987, 1988) was able to delineate the [S/Fe] trend down 
to the metallicity [Fe/H]=— 1.6 suggesting that the [S/Fe] ratio in the halo could be constant at 
a value around 0.6. The goal of this letter is to investigate whether unevolved halo stars at lower 
metallicities support the proposed trend. 

2. The Observations 

The observations were carried out in 2000 August 17-19 using the Utrecht Echelle Spectrograph 
(UES) at the Nasmyth focus of the 4.2-m WHT at the ORM (La Palma). The CCD detector 
was a 2148x2148 pixel^ SITEl. We obtained echelle spectra with a slit width ~1'.'0 which cover 
the wavelength range between 5220 and 9150 A in 20 orders at a spectral resolving power of 
R = A/AA ~ 50, 000. All spectra were background and scattered light subtracted, flat-fielded 



and extracted using standard iraf ^ procedures. The flat-field correction was attempted with a 
fast rotator, tungsten lamp and dome flats. However, the best correction from a severe fringing 
in the near-infrared was achieved with high quality dome flats. The wavelength calibration was 
performed with a Th-Ar lamp giving a dispersion of 0.087 Apixel. The final signal-to-noise ratio 
of our spectra is provided in Table 1. 

3. Spectral synthesis, stellar parameters and abundances 

The LTE spectrum synthesis code MOOG (Sneden 1973) and ATLAS9 model atmospheres 
(Kurucz 1992) were used in the present study. Atomic parameters of the sulphur doublet 8693.93 
A (log5/=-0.5, x=7.87 eV) and 8694.62 A (logc//=0.08, x=7.87 eV) were obtained from the 
VALD-2 database (Kupka et al. 1999). These logg/-values combined with the equivalent widths 
28 and 11 niA of the sulphur doublet derived from the high-resolution solar atlas (Kurucz et al. 
1984) provide a solar sulphur abundance loge(S)=7.20ib0.02 (on the customary scale in which 
loge(H)=12). This is in excellent agreement with loge(S)=7.21 given by Anders & Grevesse 
(1989). 

All our targets (Table 1) are well known metal-poor stars. The IRFM based effective 
temperature, Tofj, of two giants HD 2796, BD— 18° 5550 and a subgiant HD2665 were taken from 
Alonso, Arribas, & Martinez-Roger (1999) while the surface gravities and [Fe/H] values were 
estimated from Fe NLTE analysis (Israelian et al., in preparation) using the same tools as in 
Israelian et al. (2001). In any case, the parameters of these three evolved stars (listed in Table 1) 
are very similar to those reported in the literature (Pilachowski, Sneden &; Booth 1993, Alonso et 
al. 1999, Burris et al. 2000). The Tf.fi of five subdwarfs (Table 1) are also based on the IRFM 
and come from Israelian et al. (1998). Their metallicities and surface gravities were adopted from 
the NLTE Fe study of Thevenin & Idiart (1999) and agree perfectly with those estimated using 
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Hipparcos parallaxes (Allende Prieto et al. 1999). We note that our Tcfj-values agree with those of 
Thevenin & Idiart (1999) within 60 K. NLTE effects on sulphur are negligible (less than 0.05 dex) 
due to the weakness of the sulphur lines in our targets (Takada-Hidai et al. 2001). 

Our tests show that a typical change of +100 K in Tcfj increases [S/H] by —0.05 and —0.09 
dex for log (7=4.15 and \ogg=1.5, respectively. A change of surface gravity by +0.3 dex implies 
an uncertainty 0.11 and 0.13 dex for Tcfj=5800 and 4900 K, respectively. The errors arising from 
the uncertainties in the metallicity and the microturbulent velocity Vmic are negligible. We set 
Vmic=2 km s-^ in HD 2796, HD 2665, BD-18°5550 and Vmic='^ km s^^ for the dwarfs. The sulphur 
abundance for each star (Table 1) was obtained as the weighted mean of the abundances resulting 
for each line of the doublet (with weights proportional to the equivalent widths) . The final error 
in the sulphur abundance was derived as a quadratic sum of the errors associated with the stellar 
parameters and the errors associated to the equivalent width measurements. The final errors for 
the most metal-rich five dwarfs are similar to those reported by Frangois (1987, 1988) being of the 
order of ±0.1 dex. The errors are larger in the most metal-poor targets of the sample. 

We have reanalyzed sulphur abundances from Francois (1987, 1988) using the equivalent 
widths provided in these papers. Effective temperatures (Table 2) have been estimated using the 
IRFM calibration of Alonso, Arribas, & Martinez-Roger (1996) for Stromgren b — y and ci indices 
taken from Hauck and Mermillod (1998) and [Fe/H] adopted from Thevenin & Idiart (1999), 
Allende Prieto et al. (1999) and Allende Prieto & Lambert (1999). No reddening corrections 
were applied since all these stars are nearby within 100 parsec (Perryman et al. 1997). Typical 
errors in T^s are ±100 K. Allende Prieto et al. (1999) have shown that gravities derived from 
Fe NLTE analysis by Thevenin & Idiart (1999) agree with the gravities inferred from accurate 
Hipparcos parallaxes in the metallicity range — 1 < [Fe/H] < 0. Surface gravities of these stars 
were estimated using the Hipparcos parallaxes. We used V magnitudes from Hipparcos Catalogue 
to compute My and to estimate absolute bolometric magnitudes Mboi from Alonso, Arribas, &: 
Martinez- Roger (1995). Therefore the set of new parameters listed in Table 2 is homogeneous 
(three stars have not been considered because their parameters were rather uncertain). The masses 



of stars were inferred using theoretical isochrones from evolutionary tracks of the Geneva group 
(Schaerer et al. 1993). Our gravities have typical errors ±0.2 dex and are in good agreement with 
those presented by different authors (Allende Prieto et al. 1999, Allende Prieto & Lambert 1999) 
using, however, different evolutionary tracks. 

4. Discussion and Conclusions 

The [S/Fe] ratios obtained for our targets are listed in Table 1 and plotted in Fig. 3 together 
with Francois (1987, 1988) measurements. NLTE estimates for the metallicities of his stars have 
been made using the polynomial fit from Israelian et al. (2001) which is based on the work of 
Thevenin & Idiart (1999). Our results are consistent with those by Francois in the overlapping 
metallicity range. Remarkably, at very low metallicities we find the highest [S/Fe] ratios supporting 
a progressive increase of the [S/Fe] ratio as metallicity decreases. A minimum-square linear fit for 
the points in Fig. 3 results in a slope of —0.46 ±0.06 (l-o") which can be compared with the slope 
in the [0/Fe] vs. [Fe/H] relationship of -0.33 ±0.02 (I-ct) found by Israehan et al. (2001). 

Sulphur yields from standard Type II supernovae have been computed by many authors (e.g., 
Woosley &: Weaver 1995; Thielemann, Nomoto, & Hashimoto 1996). The [S/Fe] ratio in the ejecta 
is sensitive to the amount of mass going into the compact object formed during the supernova 
explosion (i.e depend on the adopted mass cut). There is a range of values for this mass cut which 
can perfectly provide [S/Fe] yields consistent with the observations. The tentative linear trend 
delineated in Fig. 3 may be an indication that the explosive energy of Type II supernovae may 
depend on progenitor mass. Alternatively, Nomoto et al. (2001) and Nakamura et al. (2001) 
have presented models of explosive nucleosynthesis in hypernovae. They consider exploding He 
cores with different initial masses and explosion energies in the range 1-100 x 10^^ erg. Their 
computations show that [S/Fe] ratios as those found in very low-metallicity stars can easily be 
produced in these explosions. For instance, it can be seen from Table 2 of Nakamura et al. (2001) 
that the He core with a mass 16 Mq and 10x10^^ erg explosion energy is able to provide a ratio 
[S/Fe]'~l.l in the final ejecta. It is possible that the linear trend of [S/Fe] can be explained by 
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an increasing role of hypernovae in the early epochs of the Galaxy formation (i.e., [Fe/H] < — 1). 
The estimated Galactic rate of Type Ic hypernovae is 10~^ yr^^ (Hansen 1999) but it could have 
been much higher in the early Galaxy. Small mass-loss rates due to the low metallicity in the first 
generations of massive halo stars could have led to large helium cores and the conservation of the 
primordial angular momentum. Since hypernovae are favored by stars with a large helium core 
mass and rapid rotation, a significant sulphur production in the early Galaxy could have been a 
rather common phenomenon. Similar ideas have been put forward by Qian and Wasserburg (2001) 
to explain a linear rise in the [0/Fe] ratio. Extension of their model for the sulphur case would be 
very interesting. 

Additional observational support for large amounts of sulphur synthesized during Type 
II supernova/hypernova explosion comes from the analysis of the chemical composition of the 
low-mass star orbiting the black hole in the X-ray binary system GRO J 1655-40 (Israelian et 
al. 1999). a elements are uniformly enhanced in the atmosphere of the secondary star, and in 
particular a large amount of sulphur [S/H]~0.9 was found in the secondary star. Most probably, 
this is a result of contamination from the matter ejected after supernova or hypernova explosion of 
the black hole stellar progenitor (Israelian et al. 1999, Brown et al. 2000). Very recently, Combi 
et al. (2001) have detected a bubble in the large-scale HI distribution around GRO J1655-40 
providing new evidence for a supernova explosion in this system. 

Another possible explanation for the monotonic increase in [S/Fe] at low [Fe/H] is related 
to a delayed deposition of the supernova synthesized products into the ISM due to differences in 
transport and mixing. Choosing a short mixing time (~ 1 Myr) for oxygen (a volatile), and longer 
mixing time (~ 30 Myr) for Fe (refractory), Ramaty et al. (2000, 2001) suggest that [0/Fe] should 
rise monotonically with a slope which is consistent with the observations of Israelian et al. (1998, 
2001). A similar trend is also predicted (Ramaty et al. 2001) for a volatile element like sulphur. 
As dicussed by Ramaty et al. (2000, 2001), it is not necessary that other a-elements like Mg and 
Si follow the same trend as O and S since they are refractory elements with much longer mixing 
time in the ISM. In summary, while more observations are clearly needed at very low metallicities 
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in order to confirm the trend delineated in the present work, our new measurements do not seem 
to support a [S/Fe] plateau versus [Fe/H] in halo stars. There are several scenarios that can 
account for a monotonic increase in [S/Fe] at low metallicities. 

We thank Y. Takeda and M. Takada-Hiadi for providing the NLTE corrections for sulphur 
IR doublet before the publication. This work was partially funded by project DGES PB 
98-0531-C02-02. 
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Table 1. Observing log, stellar parameters and sulphur abundances in the sample of metal-poor 

stars. The 1-a errors in the equivalent widths are computed with the formula by Cayrel (1988). 

The upper limits for HD 140283 and BD-18°555 correspond to a 3-0" values. 



Target V Exposure time (s) S/N EW (inA) EW (riiA) T^s 

8694.6 A 8693.9 A 



logs 



[Fe/H] 



[S/H] 



[S/Fe] Ref. 



HD 2665 


7.75 


3200 


400 


3.0±0.6 


- 


4990±66 


2.50±0.2 


-2.0 


1 01+0.13 
"-^•■^-^-0.14 


0.69 


1 


1 — ' 




HD 2796 


8.51 


5000 


350 


2.6±0.7 


- 


4867±58 


1.80±0.2 


-2.3 


_1 40+0.17 


0.81 


1 


HD 19445 


8.05 


10 900 


650 


2.7±0.4 


- 


5810±150 


4.46±0.2 


-1.88 


1 r,o+0.11 

-■-^•^-0.12 


0.65 


2 




HD 140283 


7.24 


2500 


300 


<2.5 


- 


5550±100 


3.74±0.2 


-2.21 


< -1.23 


<0.98 


2 




HD 157214 


5.40 


150 


350 


20±0.7 


8.0±0.7 


5625±100 


4.36±0.2 


-0.32 


-0.13±0.09 


0.19 


2 




HD 194598 


8.36 


5400 


350 


9.0±0.7 


3.0±0.7 


5960±120 


4.38±0.2 


-1.06 


-0.77±0.1 


0.29 


2 




HD 201889 


8.04 


1500 


250 


14±1 


5.0±1 


5615±100 


4.24±0.2 


-0.94 


-0.37±0.1 


0.57 


2 




BD-18°555 


9.34 


5000 


300 


<2.5 


- 


4668±63 


1.5±0.3 


-2.7 


< -1.3 


<1.4 


1 





References. — (1): Alonso ct al. (1999), Israclian ct al., in preparation, (2): Thevcnin and Idiart (1999), Israelian ct al. (1998). 
'^The signal-to-noise ratio (S/N) is provided at 8697 A. 
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Table 2. Reanalysis of sulphur abundances from the samples of Francois (1987, 1988). The 
metallicities were adopted from the references listed in the last column. Typical error in [S/Fe] is 

of the order of ±0.1 dex. 



Star Teff logg [Fe/H] [S/H] [S/Fe] Ref. 



HD 24616 


4943 


3.04 


-0.88 


-0.51 


0.37 


3 


HD 59984 


5916 


3.90 


-0.52 


-0.70 


-0.18 


1 


HD 63077 


5752 


4.17 


-0.53 


-0.55 


-0.02 


1 


HD 69897 


6263 


4.25 


-0.12 


-0.41 


-0.29 


1 


HD 76932 


5857 


4.08 


-0.79 


-0.64 


0.15 


1 


HD 88218 


5677 


3.87 


-0.42 


-0.29 


0.13 


3 


HD 91324 


6181 


3.98 


-0.23 


-0.52 


-0.29 


2 


HD 94028 


5923 


4.28 


-1.31 


-1.12 


0.19 


1 


HD 102365 


5603 


4.40 


-0.08 


-0.41 


-0.33 


2 


HD 106516 


6153 


4.30 


-0.66 


-0.56 


0.1 


2 


HD 121384 


5185 


3.39 


-0.34 


-0.79 


-0.35 


2 


HD 132475 


5619 


3.79 


-1.37 


-1.13 


0.24 


1 


HD 136352 


5637 


4.31 


-0.21 


-0.47 


-0.26 


2 


HD 148816 


5831 


4.12 


-0.51 


-0.55 


-0.04 


1 


HD 157089 


5770 


4.04 


-0.39 


-0.50 


-0.11 


1 


HD 188376 


5369 


3.62 


-0.08 


-0.30 


-0.22 


3 


HD 193901 


5668 


4.50 


-0.97 


-0.70 


0.13 


1 


HD 201891 


5841 


4.25 


-0.87 


-0.88 


-0.01 


1 


HD 203608 


6067 


4.29 


-0.57 


-0.69 


-0.12 


1 


HD 211998 


5271 


3.36 


-1.25 


-1.40 


-0.15 


1 



References. — (1): Thevenin & Idiart (1999), (2): AUende Prieto 
et al. (1999), (3): AUende Prieto & Lambert (1999). 
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Fig. 1. — Sulphur lines in the UES spectra of five metal-poor stars. 
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Fig. 2. — Synthetic spectra fits to the sulphur lines in the UES spectra (dots) of HD 19445 and 
HD 2665. The solid lines are best fits to the data corresponding to [S/Fe]=0.65 and 0.7 in HD 19445 
and HD 2665, respectively. Thin dotted lines indicate spectral synthesis with [S/Fe] shifted by ±0.15 
dex with respect the previous best fit values in each star. A vertical shift of 0.03 has been applied 
to the spectra of HD 2665 for clarity. 
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Fig. 3. — Sulphur abundances with respect to iron, derived from the near-IR S i hues. Unfilled 
squares are the data from Francois (1987, 1988) and filled circles are new measurements. The sun 
is marked with the standard symbol. Upper limits for two very metal-poor stars are not plotted. 



